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In view of the recent diboson resonant excesses reported by both ATLAS and CMS Collaborations, 
we suggest that a new weak singlet pseudo-scalar particle, ?7wzi aiay decay into two weak bosons 
while being produced in gluon fusion at the LHC. The couplings to the gauge bosons can arise 
from a Wess-Zumino-Witten anomaly term and thus we study an effective model based on the 
anomaly term as a well motivated phenomenological model. In models where the pseudo-scalar 
arises as a composite state, the coefficients of the anomalous couplings can be related to the fermion 
components of the underlying dynamics. We provide an example to test the feasibility of the idea. 


PACS numbers: 11.15.Ex, 12.60.Rc, 14.80.Va 

Introduction .— Recently, the ATLAS search for res¬ 
onant diboson {WW/ZZ/WZ) productions in fully 
hadronic final states has been reported and a discrep¬ 
ancy with the background-only model having 3 (t signifi¬ 
cance is observed around 2 TeV [1]. A similar analysis has 
been also performed by the CMS Collaboration [2] , where 
moderate excesses, less than 2cr significance, are found in 
the same mass range ~ 2TeV. In the semi-leptonic decay 
channel of the diboson resonances, however, there seems 
no excess Ian. Also, in the fully leptonic decay chan¬ 
nel of the WZ resonance, no significant deviation from 
the Standard Model (SM) prediction is observed [HIS]- 
Even if the excesses in the fully hadronic decay channels 
are not just a fluctuation, there is a possibility that a 
new resonance around 2 TeV contributes to some of the 
channels, while others are populated by misidentification 
of the boson-tagged jet: for instance, one may have a 
neutral resonance that only couples to the WW and ZZ 
channels, so that the excess in the other WZ channel is 
just a contamination owing to uncertainties of the tag¬ 
ging selections. Indeed the situation is still unclear, but it 
is worth considering the possibility for the explanation of 
the 3(7 discrepancy by the ATLAS, which is statistically 
most significant, with some new physics effects. 

This kind of exploratory exercise allows us to consider 
different new physics hypotheses and to evaluate the the¬ 
oretical motivations for various models behind such an 
experimental hint. Several authors suggest that these 
excesses may be due to the existence of a new vector 
resonance, such as a composite px or a weakly coupled 
Z' and W', and/or some other effect [THU]. Recently, 
unitarity bounds for the picture describing the excess 
with new vector resonances were also discussed 
A detailed discussion of the population of all diboson 
channels via misidentified jets can be found in Ref. [18]. 
We want to point out a novel possibility related to the 


existence of a spin-0 resonance whose couplings match 
the observed excesses while other couplings are naturally 
absent, thus giving a well-motivated phenomenological 
model. The case of a scalar coupling to gauge bosons 
via higher dimensional operators is discussed in |22| . We 
suggest instead that a new weak singlet pseudo-scalar 
particle ?7vvz: which couples to gauge bosons only via 
a Lagrangian term inspired by the Wess-Zumino-Witten 
(WZW) anomaly , can decay into two weak bosons 

[WWjZZ channels) while being produced via gluon fu¬ 
sion at the LHC. This case is theoretically well motivated 
especially in scenarios of multi-TeV scale strong dynam¬ 
ics, where such states arise as a massive scalar associated 
to an anomalous global symmetry of the confining dy¬ 
namics. In the following, we will use it as a guiding line, 
allowing the precise values of the couplings to vary in 
order to explore a larger class of models. We will rely 
on the composite scenario to establish to what extent 
the couplings can be considered “natural”, in the sense 
of the order of magnitude in the effective theory when 
particular assumptions on the specific underlying model 
are considered. We assume that this pseudo-scalar 77 .^^ 
is fermiophobic, i.e., its couplings to the SM fermions, 
in particular to the top quark, are vanishing or tiny. 
In detailed composite models this is a realistic possibil¬ 
ity and such particles with these properties are typically 
present |29j . Moreover, new scalar resonances are typ¬ 
ically expected to be lighter than their vector counter- 
partsQ so it is quite reasonable to see such a pseudo¬ 
scalar particle at a lower mass than new vectors, as in 
QCD. It is, however, required to construct an appropri¬ 
ate composite model beyond an adaptation from existing 


^ It is then easier to evade the S-parameter Eq] constraint, because 
one can push up the scale of a new strong dynamics. 
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models, because, for example, familiar dynamical models 
do not have the WW coupling [ST]. 

An obvious question is how to enhance the production 
cross section of which is typically expected to be 
tiny. Notice that the color factor = 3 counts strongly 
in TT® — >■ 77. Similarly, such an enhancement factor 
can arise from the degrees of freedom of the constituent 
fermions of which emerges via the underlying strong 
dynamics. It is therefore interesting to consider this hy¬ 
pothesis and use data from the diboson excess to have an 
idea of the coefficients required to explain such an excess, 
which in turn can hint at the more fundamental structure 
beyond the effective model. 

We therefore assume that 77-^2 couples to gluons and 
the weak bosons, as in the case of the anomaly, and does 
couple not (or only very weakly) to the SM quarks and 
leptons. Before considering detailed numbers, it is useful 
to put rough numbers on the model: the excess points 
to an effective diboson cross section of about 10 fb. This 
implies that the production cross section a-{gg —>■ r/wz) 
should be around 100 fb in order to explain the diboson 
excesses as the Branching Ratio to the WW channel is 
roughly ^ W+W-) - 2 {Nca 2 )'^/(Saj) - 10% 

when considering an anomaly induced coupling. For the 
ZZ mode, we expect a half of it. Thus the desired situa¬ 
tion, cr-Br(?7^2 W~^W~ jZZ) ^ 10 fb, can in principle 
be achieved. In this scenario, we regard the WZ excess as 
a contamination from the WWjZZ signals, noticing that 
zero events in the WZ channel maximizes a likelihood 
function in terms of the truth signal in the WWjZZjWZ 
channels m- The total width should not be so large, 
which is constrained less than, say, 100 GeV, owing to 
the one-loop effects essentially. 

One might also worry about the constraint from the 
diphoton resonance searches, because the pseudo-scalar 
can also decay into a pair of photons, unlike the vec¬ 
tor resonance. The constraints of the diphoton channel 
for a spin 0 resonance have been studied in the mass 
ranges from 150 GeV to 850 GeV by the CMS Collabo¬ 
ration [32] ■ A similar analysis was also performed by the 
ATLAS collaboration [33] • For the high-mass diphoton 
resonances, the CMS Collaborations found the constraint 
of the production cross section times branching ratio less 
than 0.3 fb for the 2TeV RS graviton |33|- The ATLAS 
Collaborations also performed a similar analysis and the 
expected -|-2cr variation limit is 0.5 fb |3S|. Even if we 
take the same bound for the spin 0 particle, our sug¬ 
gested explanation is fairly safe against this constraint in 
any case, because 77)/r(7wz ^ W+W-) - 

I2a\ ~ 0.03. The decay channel of g^z Zh is 
potentially dangerous 1311 [33, if one might expect a sim¬ 
ilar situation to the two Higgs doublet model, in which 
there appears A —>■ Zh at the tree level. In our case 
however, we can safely assume that the mixing between 
the singlet and the Higgs doublet is absent. The possible 
constraint from final states with Higgs boson(s) is thus 


easily avoided. 

In the rest of the Letter, we will present a general ef¬ 
fective description of the model, and a specihc dynamical 
model that may give rise to the desired couplings. Finally 
we will study in a model independent way the constraint 
on the couplings, necessary to reproduce the observed 
diboson excess. 

Effective Lagrangian .— The action for a weak singlet 
pseudo-scalar g^^fz with no hypercharge is 

'5'r) = y — (1) 


where is the mass of the pseudo-scalar singlet and 
the WZW term contains the effective Lagrangian for the 
diboson decay, Fwzw A / d'^x Cr^vv, witt0 
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where denotes the decay constant of the 

couplings and are arbitrary prefactors in the 

effective description, but they can be calculated in spe¬ 
cific realizations when the content of the loop terms is 
calculated, as we shall see in the next section. The cou¬ 
plings 33 , 32 and gy are, respectively, the gauge coupling 
constants of the strong, weak, and hypercharge groups. 

From the previous Lagrangian, the partial widths in 
the various channels can be easily calculated: 
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2 A term in the form Cr,WB = 
can appear through the EWSB effects, however its coefficient is 
expected to be suppressed by a I factor as it violates gauge 
invariance. 
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3 

2 

0 
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H 

3 
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n 
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2 
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d 
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2 

0 

Nl 

B 

1 

1 

0 

Nr 

B 

1 

1 
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TABLE I. Charge assignment for a vector-like model under 
the new strong dynamics SU{N), and the SM gauge sym¬ 
metries. The chirality of the fermion field is denoted by a 
subscript R or L. 


where cw = cos tw = tan 9w, e = 52 sin 9w with 9w 
being the weak mixing angle, and we define = k^+k^ 
for future reference. A naive counting of the coupling 
constants and of the numerical prefactors immediately 
shows that the production and decay to gluons will be 
the dominant channel, but the values of the couplings 
and will play a major role in the phenomenological 
results. 

This effective model allows us to easily calculate the 
diboson rates at the LHC, and check other constraints 
on the model. Before showing the numerical results, in 
the next section we will introduce a simple model of un¬ 
derlying dynamics that may lead to the required phe¬ 
nomenology. 

A Vector-like Model. — In order to discuss the expected 
phenomenology, we investigate in more detail the origin 
of the couplings and k^. In the following we 

take a simple hypothesis of a vector-like model by giv¬ 
ing the factors counting the fundamental particles in the 
anomaly loops. We do not discuss here the origin of the 
electroweak symmetry breaking, so that we assume that 
the SM-like Higgs boson with the mass being 125 GeV 
emerges as a composite object of the dynamics, or we 
incorporate it as an elementary particle. 

Let us study as an example the vector-like model shown 
in Tablejij where SU {N) represents a strongly interacting 
gauge group. Such a dynamical model with the higher 
representations of the gauge group has been studied, for 
example, in Ref. |38j . Of course, we may take the fun¬ 
damental representation as usual, if we allow arbitrary 
large N. We introduce vector-like weak doublets Q and 
L, with multiplicity nq and n^, respectively. The vector¬ 
like fermion Ni^ji is a weak singlet. The total number 
of flavors is then Nf = 2NcnQ 2nL 1, where A"c = 3 
denotes the number of ordinary QCD colors. A large 
number oi Nf is inappropriate, because the gauge the¬ 
ory loses asymptotic freedom when the fermion multi¬ 
plicity is too large. At the one-loop level, asking for 
a negative coefficient of the f3 function [3S], we obtain 
Nf < 11N/{4:T{R)), where T{R) is the trace normaliza¬ 


tion. For the two-index anti-symmetric representation, 
T{R) = (N — 2)/2. For example, the theory with N = 5 
and riQ = ul = 1 keeps asymptotic freedom. A question 
whether or not such a gauge theory might fall into the 
conformal window is beyond the perturbative approach 
and would require a dedicated analysis: some indications 
can be extracted ED]; however, only, a Lattice simulation 
[41] can give the final answer. 

The Nambu-Goldstone boson 77^2 is contained in the 
t/(l) part of the broken current SU(A^/)lX SU(Af/)ij —>■ 
SV{Nf)v. The broken current corresponding to is 

J 5 ~ - {Nf - IjiVy'^ysA^, (10) 


where we omitted the normalization factor of the axial 
current, which can be absorbed in the definition of F^j. 

We then find 

n^g=^-N{N-l)-2nQ, ( 11 ) 

i^w = \^N{N - 1) ■ {NcTiQ-\-nL), (12) 

and = 0, where A^c = 3 denotes the number 

of color. For the fundamental representation, the factor 
N{N — l)/2 should be replaced by N. The coefficient kI) 
of the WZW term for the coupling is calculated 

from the above ones and found as in this spe¬ 
cific model. The number = 2 for nq = ul = 1 

corresponds to the number of the weak doublets over that 
of the quark flavor and will play an important role in the 
diboson excess discussed later. 

For the Branching Ratios, we obtain 


Br(?7wz ^ W+W ) ^ 2{a2K^^rY 
Br(??wz 99 ) ~ 


(13) 


for nq = = 1, where we used 0:3 « 0.1 and a 2 « 0.03. 

Also, 


Bi-(??wz ^ 77) 

Br(7?wz ^ VF+1F-) 


{aKt\f 

2{a2nl)^ 



0.03, (14) 


due to in this model, where we used a = 1/128. 

These numbers can be directly compared to the exper¬ 
imental bounds on the diboson excess, and constraints 
on other channels, most notably dijet and diphoton res¬ 
onance searches: 


- ^ WW) - 10 fb, from the 

diboson excess at 2 TeV [T]; 

“ <^g9^Vwz ^ ^^{9wz 77 ) < 0-5 fb, from the 

searches of a Kaluza-Klein graviton to di-photon 
(approximate) [35] : 

- xHr(? 7^2 ^ 99 ) < 200 fb, from the search 
of dijet resonances (gluons) from a scalar [42] . 
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decay mode 

BR 

99 

83% 

WW 

11.2% 

ZZ 

3.2% 

Z7 

2% 

77 

0.4% 


TABLE II. Decay modes and branching fraction of the 
particle of 2 TeV with K^jKg = 2. 



Taking ratios of the above bounds, we can extract direct 
bounds on the ratios of Branching Ratios: 


Br(77wz ^ W+W-) 

10 

= 0.05, 

(15) 

Br(77wz ^ 99 ) 

200 

Br(77wz ^ 77 ) 

0.5 

= 0.05, 

(16) 

Br(77^2 ^ VF+1U-) 

< 10 


which are easily satisfied in this model. 

These simplified results clearly show that the fermio- 
phobic pseudo-scalar with the anomalous interactions 
can explain the diboson excesses without conflict with 
the other experimental bounds we discussed. One has to 
keep in mind, however, that a detailed model built along 
these lines may require further scrutiny concerning other 
bounds, but such a detailed study is worth pursuing only 
if the present excess will be confirmed by the ongoing 
LHC run. 

Numerical results and discussion .— In order to have 
more detailed numbers we have created a FeynRules [431 
133] model and evaluated the cross sections, branching ra¬ 
tios and decay widths numerically using Madgraph |45| . 
Using the following numerical values, ng = 1, til = 1, 
N = 2, Nc = 3, which correspond to = 2 and 
^7 ~ ~ ^?7 = 1 TeV, the production cross 

section of the 77^2 particle is 0.615 fb and its total width 
1.12 GeV at LHC with 8 TeV of center of mass energy 
for a 77^2 particle of 2 TeV of mass. 

Using instead N = 5 and all the other same parameters 
as in the previous example, increases the couplings by a 
factor of 10: = 20 and = 40, while the pro¬ 

duction cross section and width of the particle are 
a factor of 100 larger as expected (production cross sec¬ 
tion of 61.5 fb and total width of 112 GeV). The results 
for the branching fractions are given in Table |ll| These 
number are just indications based on a particular choice 
of parameters. One can see easily from the previous re¬ 
sults that increasing N (or decreasing will increase 
the cross section and allow reaching a value compatible 
with the excess. 

We consider in the following the parameters kJ in order 
to describe and bound the model in an effective way with¬ 
out reference to a particular underlying model. First, we 
can impose bounds on the couplings by taking ratios of 


FIG. 1. Cross section times branching ratios on the ACg- 
plane for Fg — 1 TeV and =0. The shaded 
region in the right upper area is excluded owing to (j{gg —> 
hwz) • Br(77wz 77) > 0.5fb. The numbers N = 4,5,6 
represent the corresponding values for the vector-like model 
with HQ = hl ~ 1. 


Branching Ratios and compare them with the bounds de¬ 
tailed in the previous section on the diboson, dijet, and 
diphoton resonant cross sections. Taking ratios of for¬ 
mulas we can eliminate the dependency on the 

cross section, and derive bounds on the couplings n]: 





1 . 45 , 


< 0.1 


92 


0.1 

sin^ 9\y 


1 . 86 , 


(17) 

(18) 


where 773 = 1.033, 772 = 0.628, and s\v? 9w = 0.2319 at 
an energy of 2 TeV. 

To compute constraints in the plane, we 

need an expression for the cross section: 


(^{99 -t ??wz) = 



(1 TeV) 2 

V 


0.615 fb 


(19) 


which can be estimated by rescaling our numerical re¬ 
sults. For the Br(77^2 99) Br(77^2 WW), by 

using Eqs. ([^-([^ for = 0, we have 


. ^ 877|(k)))2 

Br(77wz ^ 99) ^ 8gli^r + 3 g^r 


( 20 ) 


and 


Br(77wz ^ WW) 


8gl(K^,)^ + 3gl(K’^y’ 


( 21 ) 


respectively. These estimates can change if we introduce 
^WB 7^ 0 in general. In Fig. [^we show the dijet 
and diboson aww cross sections for = 1 TeV as a 
function of the and the ratio (for = 0). 

We also show the model predictions for V = 4, 5, 6. In 
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trigauge boson couplings such as for example. 

The motivation for a further more detailed analysis will 
depend on the confirmation or not of the present diboson 
excess in the near future. 
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FIG. 2. Branching ratios from the formulae of Eqs. (j^- 
Q for nPg = Q. The vector-like model with nq = = 1 

corresponds to = 2. 


the case of TV = 7, Ojj > 200 fb. The shaded region in 
Fig. corresponds to a{^gg 77 ) > 

O.Sfb. We then find that the model with TV = 4 cannot 
explain the diboson excess, on the other hand, the case 
with N = 5 does. It is fairly safe for TV = 6 , although 
the diphoton production is slightly large. The Branching 
Ratios are depicted in Fig.[^ 

In dynamical models, there usually appear many pseu- 
dos other than a singlet. How about the constraint of the 
color-octet pseudos? We estimate the difference ATVf^ of 
the mass squared by rescaling the electromagnetic mass 
splitting in the system [46] . 


- mlo \UJ a 1 ’ 


( 22 ) 


where the factor 3 is the color Casimir for the octets. We 
then find the mass of the color-octet pseudos as 3 TeV 
with Frj = 1 TeV, which is consistent with the lower 
mass bound of 2.70 TeV (2.5 TeV) by the ATLAS (CMS) 
Collaborations |32]- 

These results show that the possibility of a fermio- 
phobic pseudo-scalar singlet coupling with anomaly type 
couplings to the gauge bosons can give an explanation of 
the diboson excess without requiring the artificial sup¬ 
pression of other channels such as 77^2 The inter¬ 

pretation in terms of a more fundamental model, due to 
the large and couplings, requires relatively large 
representations as, for example, indicated in the vector¬ 
like model discussed in the previous section. This is not a 
problem in itself but the detailed model building requires 
some care in order to avoid other bounds from, for exam¬ 
ple, electroweak precision tests or the presence of other 
states which may be in the same mass range as the singlet 
77.^2 • Although there is no W^-B mixing from the model 
construction, there may appear large contributions to the 
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